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We report first evidence for a fully reconstructed decay mode of the B^r meson in the channel 
B£ -> J/^tt ± , with J/ip — > • The analysis is based on an integrated luminosity of 360 pb in 

pp collisions at 1.96 TeV center of mass energy collected by the Collider Detector at Fermilab. We 
observe 14.6 ± 4.6 signal events with a background of 7.1 ± 0.9 events, and a fit to the J/tjiir mass 
spectrum yields a Bf mass of 6285.7 ± 5.3(stat) ± 1.2(syst) MeV/c 2 . The probability of a peak of 
this magnitude occurring by random fluctuation in the search region is estimated as 0.012%. 

PACS numbers: PACS numbers: 14.40.Lb, 14.40.Nd, 13.25.Hw 



Within the standard model of elementary particles, five 
of the six different kinds of quarks combine in quark- 
antiquark pairs to make mesons. The meson com- 
bines the two heaviest of these quarks as a bottom-charm 
quark-antiquark pair. Although it has been observed in 
semileptonic decay modes [l|, |2| , up to now no evidence 



for the B^ has been found in fully reconstructed decay 
modes |j. Lit [[J |(| Consequently, its mass M(B C ) has not 
been measured with good precision. 

Nonrelativistic potential models predict the b and c 
quarks to be tightly bound with a ground state mass 
in the approximate range 6200-6300 MeV/c 2 0, S 0. 
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Recent QCD-based perturbative computations up to 
0{a 4 s ) predict M(B C ) to be 6307 ± 17 MeV/c 2 H3, EH 
Most recently, a three-flavor lattice QCD calculation ob- 
tains M(B C ) = 6304 ± 12(stat © syst) ±J 8 (cutoff ef- 
fects) MeV/c 2 0. 

Several of the predicted Bf decay modes contain a 
J/ip meson 0]. These are among the most easily re- 
constructable Bf decays at CDF, owing to an efficient 
dimuon trigger giving high purity J/ip — > fi + n~ recon- 
struction. The CDF collaboration made the first observa- 
tion of the Bf meson in the semileptonic decay channels 
Bf — ► J/ipl^UiX, in a sample of 110 pb _1 of data at 
y/s = 1.8 TeV in Run I at the Tevatron Q. The symbol 
X denotes possible undetected decay particles. With a 
signal of 20.4+gJ events, the Bf mass was measured to 
be 6.40 ± 0.39(stat) ± 0.13(syst) GeV/c 2 . Recently, the 
DO Collaboration reported a preliminary observation of 
a Bf signal in the decay channel B^_—> J / 'ip ^ v in a 
sample of 210 pb _1 of Run II data [2j. 

In this paper we report first evidence for the Bf meson 
in the fully reconstructed decay channel Bf — ► J/ipn^, 
with J/ip — ► The analysis is based on a dataset 

of 360 pb _1 in pp collisions collected at y/s — 1.96 TeV 
by CDF at the Tevatron during Run II. 

The CDF II detector consists of a magnetic spectro- 
meter surrounded by calorimeters and muon chambers 
and is described in detail elsewhere (l4|. The compo- 
nents most relevant to this analysis are briefly described 
here. The tracking system is in a 1.4 T axial magnetic 
field and consists of a silicon microstrip detector (LOO, 
SVX, ISL, in increasing order of radius) 0,0,0] sur- 
rounded by an open-cell wire drift chamber (COT) 0. 
The muon detectors used for this analysis are the central 
muon drift chambers (CMU), covering the pseudorapid- 
ity range \r)\ < 0.6 [ilj. Eof . and the extension muon drift 
chambers (CMX), covering 0.6 < \rj\ < 1.0. Cylindrical 
coordinates are used with the +z axis in the proton beam 
direction. 

This measurement uses events containing pairs of 
muons, each with | ^7 1 < 1.0, selected with a three- level 
trigger. At the first trigger level, muon-candidate track 
segments in CMU and CMX are matched to COT tracks 
obtained with a hardware processor |2l| . Dimuon trig- 
gers use combinations of CMU-CMU and CMU-CMX 
muons with p T >1.5 (2.0) GeV/c for CMU (CMX) 
muons, where px is the momentum transverse to the 
beamline. At the second level, opening angle and 
opposite-charge cuts are imposed on the muon pairs. 
At the third level, three dimensional (3-D) tracking is 
performed to select muon pairs with invariant mass, 
M{n + iT), between 2700 and 4000 MeV/c 2 . 

To reconstruct the Bf — > J/ipn^ decay offline, we 
make several requirements on the quality of the tracks 
and the J/ip candidate. To ensure good primary and sec- 
ondary vertex resolution, each track must have an r — <p 
position measurement on at least three of five SVX layers. 



For J/ip identification, we require matching between the 
COT muon tracks and the muon chamber track segments. 
In addition, we require that 3042 < M(/i + /i~) < 3152 
MeV/c 2 , the average J/ip mass resolution in our sam- 
ple being 14 MeV/c 2 . Each other charged particle track 
with pt > 400 MeV/c is treated as a pion candidate to 
be combined with the J/ip. The pion candidate and the 
two muons are then fitted to a common 3-D vertex, with 
M(u + U~) constrained to the world average J/ip mass 
value [22| . All combinations for which the vertex fit con- 
verged are retained. The primary vertex position is cal- 
culated from the other tracks in each event. 

The Bf search was performed using the following anal- 
ysis method. The mass values of the J/ipir combina- 
tions in the search window 5600 < M(J/-07r ± ) < 7200 
MeV/c 2 , referred to as Bf candidates, were temporarily 
hidden. The search window was chosen to correspond to 
the ±2 standard deviation region around the CDF Run I 
measurement of the Bf mass 0; it is approximately 100 
times wider than the expected Bf mass resolution. 

In order to optimize the significance of a possible Bf 
signal, we varied the selection criteria to maximize the 
function Q — Sf/0-5 + y/B av ) [2j|. Here, Sf is the ac- 
cepted fraction of signal events, in this case taken from 
a Monte Carlo sample, and the background B av is the 
number of selected Bf candidates within the search win- 
dow, scaled to correspond to a mass range of 63 MeV/c 2 , 
based on the average mass resolution of a Bf candidate 
within the search window. The term 1.5 is appropriate 
for optimizing a search for a signal at least 3cr above 
background fluctuations. The distributions of the selec- 
tion variables for the signal events were evaluated using 
samples of simulated Bf — > J/ipn^ decays. These were 
generated with a Bf mass of 6400 MeV/c 2 , a lifetime of 
0.46 ps 0, and pr and rapidity distributions according 
to a leading order perturbative QCD calculation [24|. A 
harder px spectrum 0] was used as an alternative to 
check the stability of the optimal selection criteria; these 
were not very sensitive to variations of the px spectrum 
or the assumed lifetime within its experimental uncer- 
tainty. The Monte Carlo Bf decays were processed with 
full detector simulation and the same trigger and recon- 
struction criteria as the data. The distributions of the 
selection variables for the background were taken from 
the data in the search window, in which the contribution 
from a signal is expected to be small. 

Optimized cuts were determined for the following se- 
lection variables: the J/tpn^ 1 three-track 3-D vertex fit 
(% 2 < 9 for four degrees of freedom), the pion track con- 
tribution to the vertex fit (% 2 < 2.6), the impact param- 
eter in r — cp of the Bf candidate with respect to the 
primary vertex (< 65 /^m), the maximum ct where t is 
the proper decay time of the Bf candidate (< 750 /^m), 
the transverse momentum of the pion (> 1.8 GeV/c), the 
3-D angle between the momentum of the Bf candidate 
and the vector joining the primary to the secondary ver- 
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FIG. 1: The invariant mass distribution of the B — > J/i/jK' 
candidates. The curve is a fit to the data. 



tex ((3 < 0.4 rad), and the significance of the projected 
decay length of the Bf candidate onto its transverse mo- 
mentum direction (L xy /a(L xy ) > 4.4). After these selec- 
tions, 390 candidates remain in the search window, with 
no two candidates from the same beam crossing. 

A sample of i? ± mesons, reconstructed in the decay 
mode B^ — > J/ipK^, was analyzed as a control sam- 
ple in order to check our understanding of the recon- 
struction of the relevant variables in the simulation. The 
B^ — > J/ij) K ± decay topology is the same as that of 
B^r — > J/ip tt^~ , apart from the different masses and 
lifetimes. The B mass distribution, shown in Fig. 1, 
was obtained using the same selection requirements as 
optimized for the Bjf candidates, but without the cut 
on maximum ct. A total of 2378 ± 57 B ± -> J/V>if ± 
signal events is found, with a fitted mass of 5279.0 ± 0.3 
MeV/c 2 . The fit takes into account a small contribution 
from the Cabibbo-suppressed decay S ± — > J/V'7r ± . The 
average mass resolution is 11.5 ± 0.3 MeV/c 2 , in agree- 
ment with the simulation, which can thus be used with 
confidence to evaluate the expected mass resolution for 
Bf decays. The B ± yield is used to calculate the ex- 
pected B^r yield. The relative trigger and reconstruction 
efficiency, e B ±/es±, is in the range 35%-85%, with un- 
certainties arising from the Bf- px spectrum and the Bf 
lifetime. On the basis of the B ± yield, previous CDF 
cross section measurements llj, and theoretical calcula- 
tions [IlliimiHlilliSliJof the branchin g fractions 
of the B^r — > J/tpir^ and B^r — > J/tpl u decay modes, a 
B^ yield in the range of 10 to 50 events is expected. 

A search procedure was then defined to identify any 
possible signal in the data and to estimate its signifi- 
cance. This was based on a scan of the search region in 
10 MeV/c 2 intervals, with a sliding fit window extend- 
ing from -100 MeV/c 2 to +200 MeV/c 2 in mass around 



each nominal peak position, m. This window was chosen 
to minimize possible contributions from partially recon- 
structed Bf decays below the peak position (e.g. into 
J/ip and more than one additional particle). For each 
value of m, a fit function was defined as a Gaussian sig- 
nal with mean m, combined with a linear background 
term. The Gaussian resolution was fixed as a linear func- 
tion of m based on Monte Carlo simulation, and varied 
from 13 to 19 MeV/c 2 over the search region. The three 
fit parameters were the number of signal (S) and back- 
ground (B) events and the linear background slope. The 
output of a scan was defined to be the largest value of 
£ = 5/(1. 5 + \/B), S maa; , obtained from the 131 fits per- 
formed in the mass interval 5700 < M{ J/ip^) < 7000 
MeV/c 2 . 

The distribution of T, max for the null hypothesis was 
obtained from Monte Carlo experiments [32], in which 
the mass spectra were derived from a smooth background 
model. This model was necessarily approximate owing to 
the initially hidden mass distribution. The model con- 
sisted of a linear background, to describe combinatoric 
events, and a "physical" background to describe partially 
reconstructed B^r decays in the mass range below 6400 
MeV/c 2 . Studies showed that the main source of combi- 
natoric background are events in which a genuine J/psi 
is paired with an uncorrelated track. The shape of the 
physical background was based on Monte Carlo simula- 
tions of inclusive Bf — ► J/ipX decays, with branching 
ratios taken from Ref. |l3j ]. 

Applied to the 390 event data sample, the scan pro- 
cedure found a T, max near m = 6290 MeV/c 2 , which is 
compatible with a B^r signal of 19 ± 6 events. Using 
a large set of Monte Carlo simulations we modeled the 
shape of the observed background, and, analysing it in 
the same way as the data, evaluated the probability that 
a random enhancement has a T, max value exceeding that 
of the data. This probability was found to be 0.17%. 

After the above steps had been performed, further 
checks on the previously hidden events revealed that the 
existing pion selection allowed two classes of fitted tracks 
that were unsuitable for the Bf search. The first class 
had insufficient number of COT hits to give good mass 
resolution and so was not compatible with a search for a 
narrow Gaussian signal; the second class had poor SVX 
resolution in the z direction and was dominated by com- 
binatorial background. The above two classes of events 
contributed 10% to the signal; they would be ex- 
pected to contribute fewer than two events to the B^r 
signal, but they increase the combinatorial background 
by about 40% over the J/tpn^ mass range. After re- 
moval of both classes of poor quality tracks in addition 
to the original optimized cut selection, 220 candidates re- 
mained. These were required to have good SVX z resolu- 
tion on both the pion track and at least one of the muon 
tracks. This final track selection, which maximizes Q, is 
therefore not fully blind; it is based also on the observed 
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properties of the B ± signal and the overall properties of 
the B~- candidate sample. 

Figure 2 shows the mass spectrum for the 220 event 
sample. The main features are the B^r — > J/tpn^ sig- 
nal peak near 6290 MeV/c 2 , a linear combinatorial back- 
ground above this peak, and a broad enhancement below 
the peak which can be attributed to the physical back- 
ground from partially reconstructed B^ decays. We per- 
form a global unbinned likelihood fit over the entire mass 
range to obtain the mass and yield for the Bf signal. The 
fit included a Gaussian signal with a variable mass but 
with a resolution whose mass-dependent value was de- 
termined by the Monte Carlo simulation, together with 
background modeled as a linear combinatorial term and 
a broad low-mass Gaussian contribution for the physical 
background. A signal of 14.6 ± 4.6 events is obtained 
centered at a mass of 6285.7 ±5.3 MeV/c 2 . The stan- 
dard deviation of the Gaussian at the central value of the 
signal mass is 15.5 MeV/c 2 . The background within a re- 
gion of ±2 standard deviations from this mass value is 
7.1 ± 0.9 events. The statistical significance of the signal 
is discussed below. Within the signal region, the distri- 
butions of the selection variables agree within statistics 
with those of the Monte Carlo simulation. 

Systematic uncertainties on the B^r mass determina- 
tion due to measurement uncertainties on the track pa- 
rameters (±0.3 MeV/c 2 ) and the momentum scale (±0.6 
MeV/c 2 ) are evaluated from the corresponding uncer- 
tainties on the B^ mass analysis [3^ | . Further uncertain- 
ties are due to the possible differences in the pt spectra 
of the B^ and Bf- mesons (±0.5 MeV/c 2 ) and our lim- 
ited knowledge of the background shape used in the final 
mass fit as well as uncertainty in the signal width (±0.9 
MeV/c 2 ) [34j. The total systematic uncertainty is evalu- 
ated to be ±1.2 MeV/c 2 . 

The signal peak is robust under variations of the pion 
track quality selection. We have investigated several 
methods for determining the best figure of significance 
for such a peak over a broad mass range. The method 
that gives the best sensitivity to a real signal is based on 
the standard significance measure S/VB. We repeated 
the Monte Carlo scans for the new track selection to de- 
termine the null hypothesis distribution for S/^fB. Ap- 
plying to the Monte Carlo simulations the same global fit 
method as to the data, we find that the probability that 
a random enhancement anywhere in the range 5800-7000 
MeV/c 2 exceeds the value of S/^/B for the experimental 
peak is 0.012%. 

In view of the limited statistics of the observed 
mass peak, an independent consistency check was per- 
formed. If the mass peak is due to fully reconstructed 
Bf — > J/ip^ decays, partially reconstructed Bf — > 
J/V'+track+A decays should be detectable in the mass 
region below the peak but not in the region above. The 
pion candidate in partially reconstructed decays should 
have a small impact parameter d xy relative to the J/tp 
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FIG. 2: The invariant mass distribution of the J/ip^ can- 
didates and results of an unbinned likelihood fit in the search 
window. The inset shows the peak section of the distribution. 
The broad enhancement below 6.2 GeV/c 2 is attributable to 



partially reconstructed B^r mesons. 



vertex, consistent with being physically associated with 
it, whereas the pion candidate in combinatorial back- 
ground events should have a broad d xy distribution re- 
flecting random association with the J /ip vertex. 

To investigate this, we relax the cuts on (3, the impact 
parameter of the Bf candidate, and the \ 2 of the 3-D 
vertex fit, so as to make a signal in the d xy distribution 
visible over the broader combinatorial background. We 
compare the distribution of d xy of the pion candidate 
in the region 5600 < M(B C ) < 6190 MeV/c 2 (lower side 
band) to that in the region 6390 < M(B C ) < 7200 MeV/c 2 
(upper side band), where the main contribution should 
be combinatorial. 

Figure 3 (top) shows the difference between the lower 
(4900-5100 MeV/c 2 ) and upper (5400-5700 MeV/c 2 ) 
sidebands for the d xy distribution in the B^ data sam- 
ple, with a large excess of events visible at small d xy 
values. Figure 3 (bottom) shows the corresponding plot 
obtained using the B^r candidate sample. An enhance- 
ment is visible with a shape compatible with that seen 
in the S ± sample. The B ± curve, rescaled to fit the 
B^r data, provides a good description of this distribu- 
tion. The excess of low d xy events in the B^r sample is 
evaluated to be 244 ± 59, where the uncertainty is sta- 
tistical only. This result is consistent with Monte Carlo 
estimates based on the calculations of 13] . This supports 
the hypothesis that the broad physical background below 
the signal peak, evident in Fig. 2, is in fact associated 
with partially reconstructed Bf decays. 

In conclusion, we observe a peak in the J/tpir mass 
spectrum at a mass of 6285.7 ± 5.3(stat) ± 1.2(syst) 
MeV/c 2 . This peak is consistent with a narrow, weakly 
decaying particle state and is interpreted as the first ev- 
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FIG. 3: Impact parameter of the third track relative to the 
J/'ijj vertex for the lower sideband region, after subtraction of 
the same distribution for the upper side band: (top) the curve 
is the sum of two Gaussians, fitted to the B data points; 
(bottom) the B^r data points, overlaid with the above curve, 
rescaled. In both cases the selection criteria were relaxed. 



idence for fully reconstructed decays of the meson. 
The mass value has much improved precision over the re- 
sults obtained in B^r semileptonic decays P, Q . There is 
also good agreement with recent theoretical predictions 
for the Bf mass around 6300 MeV/c 2 HUE], El 

We thank the Fermilab staff and the technical staffs of 
the participating institutions for their vital contributions. 
We also thank A.V. Berezhnoy, C.H. Chang and X.G. Wu 
for making available their calculations of B^ production 
spectra. This work was supported by the U.S. Depart- 
ment of Energy and National Science Foundation, the 
Italian Istituto Nazionale di Fisica Nucleare, the Ministry 
of Education, Culture, Sports, Science and Technology of 
Japan, the Natural Sciences and Engineering Research 
Council of Canada, the National Science Council of the 
Republic of China, the Swiss National Science Founda- 
tion, the A. P. Sloan Foundation, the Bundesministerium 
fiir Bildung und Forschung, Germany, the Korean Sci- 
ence and Engineering Foundation and the Korean Re- 
search Foundation, the Particle Physics and Astronomy 
Research Council and the Royal Society, UK, the Russian 
Foundation for Basic Research, the Comision Intcrminis- 
terial de Ciencia y Tecnologia, Spain, in part by the Eu- 
ropean Community's Human Potential Programme un- 
der contract HPRN-CT-2002-00292, and the Academy 
of Finland. 



[2 

[3 

[4 

[5 

[6 

[7; 
[8 
[9 
[10 

[11 
[12 

[13. 
[14 

[IS 
[16 

[17 

[is; 

[19 

[20 

[21 
[22' 
[23 
[24 

[25' 



[26 
[27 

[28 

[29 
[30 
[31 
[32 



[33 
[34 



(1998) . 

D. Lucchesi, in ICHEP 2004, ed. H. Chen et al, World 
Scientific, Singapore, 2005, pl58. 

CDF Collaboration, F. Abe et al, Phys. Rev. Lett. 77, 
5176 (1996). 

DELPHI Collaboration, P. Abreu al, Phys. Lett. B398, 
207 (1997). 

ALEPH Collaboration, R. Barate et al, Phys. Lett. 
B402, 213 (1997). 

OPAL Collaboration, K. Ackerstaff et al, Phys. Lett. 
B420, 157 (1998). 

W. Kwong and J. Rosner, Phys. Rev. D44, 212 (1991). 

E. Eichten and C. Quigg, Phys. Rev. D49, 5845 (1994). 
S. Godfrey, Phys. Rev. D70, 054017 (2004). 

N. Brambilla, Y. Sumino and A. Vairo, Phys. Rev. D65, 

034001 (2002). ; 

N. Brambilla et al, hep-ph/0412158 and references 
therein. 

I.F. Allison et al, Phys. Rev. Lett. 94, 172001 (2005) 
and Nucl. Phys. Proc. Suppl. 140:440-442 (2005). 
V.V. Kiselev, Phys. Atom. Nucl. 67, 1559 (2004). 
CDF Collaboration, D. Acosta et al, Phys. Rev. D71, 
032001 (2005). 

A. Sill et al, Nucl. Instrum. Methods A447, 1 (2000). 
A. Affolder et al, Nucl. Instrum. Methods A453, 84 

(2000) . 

C.S. Hill et al, Nucl. Instrum. Methods A530, 1 (2004). 
T. Affolder et al, Nucl. Instrum. Methods A526, 249 
(2004). 

G. Ascoli et al, Nucl. Instrum. Methods A268, 33 
(1988). 

T. Dorigo et al, Nucl. Instrum. Methods A461, 560 

(2001) . 

E.J. Thomson et al, IEEE Trans. Nucl. Sci. 49, 1063 

(2002) . 

S. Eidelman et al, Particle Data Group, Phys. Lett. 
B592, 1 (2004). 

G. Punzi, Proceedings of PhyStat2003, (SLAC, 2003), 
eConf C030908, MODT002 (2003). 

C.H. Chang and X.G. Wu (private communication); C.H. 
Chang et al, CPC, 159 (2004) 192; C.H. Chang and X.G. 
Wu, Eur. Phys. J. C. 38 (2004) 267. 

A.V. Berezhnoy (private communication); A.V. Berezh- 
noy, V.V. Kiselev, A.K. Likhoded, Z.Phys.A356, 79 

(1996) ; A.V. Berezhnoy, A.K. Likhoded, M.V. Shevlya- 
gin, Phys.Atom.Nucl.58, 672 (1995). 

C.H. Chang et al, Phys. Rev. D49, 3399 (1994). 

A.V. Berezhnoy, et al, Phys. Atom. Nucl. 60, 1729 

(1997) . 

A.Yu. Anisimov, et al, Phys. Atom. Nucl. 62, 1739 

(1999) . 

P. Colangelo et al, Phys. Rev. D61, 034012 (2000). 
A. A. El-Hady et al, Phys. Rev. D62, 014019 (2000). 
K. Anikeev et al., |hep-ph/0201071| 

WA. Rolke and A.M. Lopez, Proceedings of Phy- 
Stat2003, (SLAC, 2003), eConf C030908, MODT002 

(2003) . 

CDF Collaboration, D. Acosta et al, |hep-e x/0508022 
submitted to Phys. Rev. Letters. 

L. Nicolas, Ph.D. dissertation, Univ. of Glasgow, 2005. 



[1] CDF Collaboration, F. Abe et al, Phys. Rev. Lett. 81, 
2432 (1998); F. Abe et al, Phys. Rev. D58, 112004 



